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We show that a chiral / + i/-wave superconducting pairing may be induced in the lowest heavy 
hole band of a hole-doped semiconductor thin film through proximity contact with an s-wave super- 
conductor. The chirality of the pairing originates from the 3ir Berry phase accumulated for a heavy 
hole moving along a close path on the Fermi surface. There exist three chiral gapless Majorana edge 
states, in consistence with the chiral / + z/-wave pairing. We show the existence of zero energy 
Majorana fermions in vortices in the semiconductor-superconductor heterostructure by solving the 
Bogoliubov-de-Gennes equations numerically as well as analytically in the strong confinement limit. 
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Superconductors/superfluids with unconventional 
pairings have been an important subject in condensed 
matter physics for many decades because of their rich 
physics and important applications. There has been 
considerable experimental evidence to support that the 
pairing symmetry in high-T c superconductors is the 
d-w&ve The pairing symmetry in the superfluid 

3 He systems has been found to be p-wave Q- The 
superconducting order parameters in Sr2Ru04 and 
some heavy-fermion materials are suggested to be the 
chiral p x + ip y wave Q, although the true nature of 
the order parameters in these materials has not been 
fully settled in experiments [1, The importance of 
the chiral p-wave superconductor/superfluid is that the 
quasiparticle excitation inside a vortex is a zero-energy 
Majorana fermion with non-Abelian exchange statistics, 
which is a crucial ingredient for topological quantum 
computation (TQC) [6j. 

However, in contrast to the simple s-wave supercon- 
ductor described by the BCS theory, the theoretical 
description and experimental identification of uncon- 
ventional superconductivity and the associated exotic 
physics in natural solid state systems are often difficult, 
and in many systems, controversial. For instance, despite 
the tremendous technological potential, the observation 
of the exotic properties such as quantum half-vortices 
and non-Abelian statistics in Sr2Ru04 has been a seri- 
ous problem because of the small quasiparticle excitation 
energy gap as well as the intrinsic spin-orbit coupling in 
the suggested p-wave order parameter Q. Therefore it 
should be not only interesting, but also important to in- 
vestigate whether various unconventional superconduct- 
ing pairings and the associated exotic physics can be ex- 
ternally induced from conventional s-wave superconduc- 
tors/superfluids [sl— Tlllj . For instance, schemes have been 
proposed recently to induce the zero-energy Majorana 
fermions in the vortex core of conventional s-wave super- 
conductors that are proximately coupled to topological 
insulators or electron-doped semiconductors 0- 11 . 



In this paper, we propose that a chiral / + i/-wave 
superconducting pairing can be induced in a hole-doped 
semiconductor thin film through the proximity contact 
with an s-wave superconductor. To the best of our knowl- 
edge, a chiral / + i/-wave superconducting pairing and 
the associated exotic physics have not been unequivocally 
identified in any condensed matter system. The induced 
chiral / + i/-wave pairing symmetry has a topological 
origin: the geometric phase 12 , l]| of holes in the Bloch 
band. It is well known that an electron/hole evolving 
adiabatically in the reciprocal space accumulates a geo- 
metric phase associating with the adiabatic change of the 
quasi- momentum (l3j , in analogy to the Aharanov-Bohm 
phase acquired by electron moving in the real space in 
the presence of a magnetic field. The geometric phase 
is nonzero in the hole-doped semiconductors with non- 
vanishing spin-orbit coupling, which tunes an original s- 
wave pairing into a chiral / + i/-wave pairing for holes in 
the lowest energy band. The induced chiral / + i/-wave 
superconductor has a full pairing gap in the 2D bulk, 
and 3 gapless chiral Majorana fermions at the edge. By 
solving the Bogoliubov-de-Gennes (BdG) equations an- 
alytically and numerically, we show that there exists a 
Majorana zero energy state in the vortex core of the 
semiconductor-superconductor heterostructure in some 
parameter regions. The corresponding quasiparticle ex- 
change statistics in this system is the same as that for a 
chiral p-wave superconductor or superfluid, therefore the 
proposed heterostructure can be used as a platform for 
observing non-Abelian statistics and performing TQC. 
The advantage of using hole-doped, instead of electron- 
doped, semiconductors for TQC is that hole-doped semi- 
conductors have stronger spin-orbit coupling due to the 
larger effective mass of holes and the p-like symmetry of 
the valence band, resulting in larger carrier densities. 

The physical system we consider is a heterostructure 
composed of an s-wave superconductor, a hole-doped 
semiconductor thin film, and a magnetic insulator (Fig. 
[1^). In the semiconductor thin film, the dynamics of 
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Figure 1. (Color online) (a) A schematic illustration of the 
heterostructure composed of a hole-doped semiconductor thin 
film, an s-wave superconductor, and a magnetic insulator, (b) 
The Band structure of the hole-doped semiconductor. The 
lower (upper) two bands are the HH (LH), respectively, a — 
2 x 10 5 m/s, h-s/Jj^j = 3 x 1(T 26 kg-m/s, 71 = 6.92, 72 = 2.1, 
fto = 1-75 meV. The parameters are chosen from GaAs (other 
materials such as InAs, InSb yield the same physics). The 
hole density is 



holes can be described by a single particle effective Hamil- 
tonian that contains both Luttinger four band model and 
spin-3/2 Rashba term [l4| 



Ha = 



(71 + 572/2) k 2 - 2 72 (k • J) 2 J /2m 
+a (J x k) • z + 2h J z - fi, (1) 



where J is the total angular momentum operator for a 
spin-3/2 hole, 71 and 72 are the Luttinger parameters, 
/1 is the chemical potential. Henceforth, we have set 
H = 1. The confinement of the quantum well along the 
z direction makes the momentum be quantized on this 
axis, that is, (k z ) m 0, (fc 2 ) = (tt /a) 2 , where a is the 
thickness of the quantum well, a is the Rashba spin- 
orbit coupling strength. The crucial difference between 
the Rashba terms in the 2D hole and electron gases is 
that J in the 2D hole gas is a spin-3/2 matrix, describ- 
ing both the heavy holes (HH) and light holes (LH) . The 
term 2HqJ z describes a Zeeman splitting induced either 
through the polarization of the local magnetic moments 
in the semiconductor [l6| or the exchange field through 
the contact with a magnetic insulator. 

The eigenstates of the Hamiltonian ([lj can be written 

as 



<T/k = ( u e 
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where #k is the azimuthal angle of k, Ui are functions of 
k only and the eigenstates of the reduced Hamiltonian 

H = - l2 k 2 J 2 y /m - akJ x - l2 (k 2 z ) J 2 /m + 2h J z . (3) 

The particular choice of the wavefunction in ([2]) ensures 
that the wavefunction for the lowest HH band is single- 
valued at k = (i.e., only 113 ^ at k = 0). Additional 
total phase e" n6k need be multiplied to the wavefunctions 
for the other bands to ensure the single value. In Fig. 
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Figure 2. (Color online) (a) Plot of A(k) with respect to 
k. (b) Plot of the bulk quasiparticle energy E n of the BdG 
equation ([9]). Solid, dashed, dotted, and dashed dotted lines 
correspond to the lowest to the highest bands in Fig. [TJj. 
fi — —32.5 eV. The other parameters are the same as that in 
Fig. 1. 



Qi>, we plot the eigenenergies with respect to k for the 
2D hole gas. The degeneracy between different HH and 
LH bands at k = is lifted due to nonzero (fc 2 ) and the 
Zeeman field ho. In the strong (weak) confinement region 
272 (k 2 ) I'm > (<) 4/io, the second lowest energy band is 
the HH (LH), and the corresponding wavefunction of the 
band can be written as ^ k e' l3ek (^ke I0k ). 

The proximity-induced superconductivity in the hole- 
doped semiconductor can be described by the Hamilto- 



17] 



H„ 



E 



mj=l/2,3/2 



J dr{A s (r) cl ]C l mj +H.c] 



(4) 

where cj^ are the creation operators for holes with the 
angular momentum mj and A s (r) is the proximity- 
induced gap. When the chemical potential lies between 
the lowest two bands (Fig. [TJd), and only the lowest HH 
band is occupied with holes, the effective superconduct- 
ing pairing for holes becomes 



(a-kd-k) oc iA s g (k) exp (i'S9u) 



(5) 



where Ok = X)m, Xm./Ckm., is the annihilation operator 
for holes at the lowest HH state ^k, with the coefficient 
Xmj = *kmj) 9{k) = 2(u u 3 + uiu 2 ). We have used 

A s oc (akmjfl-k(-mj)) and (akmjfl-kra'j) = if mj ^ 
—m'j to derive Eq. ©. Clearly the pairing order A e // 
has a chiral / + if -wave symmetry. Around the Fermi 
surface g (k) — > 1, that is, A e // — > ?A S exp (i3#k)- 

The 30k phase in A e ff originates from a 37r Berry 
phase accumulated when the holes move in the momen- 
tum space. In the lowest HH band, the Berry phase for 
a hole along a loop on the Fermi surface is 



A-dk = A{k F )59 k , 



(6) 



ki 



where the Berry connection A = i ("Fkl Vk l^k) 
A (k) V0 k with A (fc) = (3u 2 + 2u\ + uf) , 66^ = ka — 



is the change of the azimuthal angle from ki to k2. In 
Fig. [2^, we see A(k) — > 3/2 around the Fermi surface 
for the lowest HH state, indicating a 359 k /2 Berry phase 
(37r for a close loop) for a single hole and 359^ phase 
for a Cooper pair. Therefore A e // in the lowest band 
has a phase factor exp (z3#k)- Similarly, we can calculate 
the Berry connection for the other three upper bands. 
In Fig. [2^, we see the corresponding Berry phases are 
— 37r, 7r, and —tt, respectively, which means that the su- 
perconducting pairing symmetries for holes at the upper 
HH, lower LH, and upper LH bands are chiral / — if, 
Px + ip y , Px — ipy-w&ves, respectively. 

The physics origin of the chiral / + if -wave supercon- 
ducting pairing in the lowest HH band is more transpar- 
ent in the strong confinement limit (fc <C V\fcf}), where 
the four band Hamiltonian ([T]) can be diagonalized into 
two effective two-band Hamiltonians for the HH and LH 
respectively. The effective Hamiltonian Hhh for the HH 

Hhh = ri k A + rnk 2 + il3 (k 3 _a + - k\aJ) +3h Q a z - p.. (7) 

Here k± = k x ± ik y , a± = (a x ± ia y ) /2 are Pauli ma- 
trices applied on the two HH states (denoted as pseu- 
dospin | and i), r\i are the reduced coefficients, /3 is 
the effective coupling strength, /2 is the effective chem- 
ical potential. The Hamiltonian ([7]) is similar as the 
Rashba type of Hamiltonian for electron-doped semi- 
conductors except that k x + ik y is now replaced with 
(k x + iky) 3 = A: exp (i3$k) and there is a fc 4 term to en- 
sure the bands bend up for a large k. Therefore a chiral 
/ + if-w&ve superconducting pairing is obtained when 
only the lowest HH band is occupied with holes. 

A chiral / + i/-wave superconductor should have a full 
pairing gap in the 2D bulk, and C = 3 gapless chiral 
Majorana fermions at the edge (l8j |. where 



2n ^-E„<o 



d 2 k m 



(8) 



is the first Chern number in the momentum space, fi™ = 
—2Im(d& n /dk x \d& n /dky) is the Berry curvature of the 
n-th band. E n and <!>„ are eigenenergies and wave func- 
tions of the BdG equation, which, in the Nambu spinor 
basis, can be written as 



H 



Mr) 



A*(r) -a x T y H^T y a x 



$„(r) - EA(r). (9) 



tion, a x = diag(a x ,a x ), r y 



In a 



Here $„(r) = [u„,3/2 U n ,l/2 Un,-l/2 U n ,-3/2 «n,-3/2 
v n.-i/2 ~ v n,i/2 ~ v n,3/2] T is the quasiparticle wavefunc- 

-ihx2 
ihx2 

uniform system with a constant A s (r), the BdG equa- 
tion (|9|) can be solved in the momentum space and the 
corresponding quasiparticle energy dispersions E n (fc) are 
plotted in Fig. [5b. We see a 2A S energy gap is opened 
at the Fermi surface. Using the eigenwavefunctions $„ 
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Figure 3. (Color online) Plot of the wavefunction of the zero 
energy state, fi = —32.5 meV lies in the gap between two HH 
bands. The other parameters are the same as that in Fig. 1. 
Inset: Plot of the s-wave pairing gap with a vortex. 



in the lowest four bands of the BdG equation (|9]) (i.e., 
E n < 0), we confirm that the Chern number C = 3, 
which is consistent with the chiral / + if -wave super- 
conducting pairing and yields 3 gapless chiral Majorana 
fermions at the edge of the superconductor. 

The chiral / + if -wave pairing may lead to novel ex- 
otic physics that has not been explored before (e.g., 
fractional Josephson effects [H|]). Here we focus on 
the Majorana fermions in vortices in the heterostruc- 
ture that can be used for TQC. In the presence of a 
vortex in the heterostructure, the pairing order param- 
eter takes the form A s (r) = A s (r)e l9 . For simplic- 
ity of the calculation, we consider a 2D cylinder ge- 
ometry with a hard wall at the radius r — R and 
a single vortex at r — 0. This system preserves the 
rotation symmetry and the BdG equation can be de- 
coupled into different angular momentum channels in- 
dexed by I with the corresponding spinor wavefunction 
<(r) = e afl [< a/2 e- w 



p iS ,.l 
,-l/2 c ' u n,-3/2 L ' 



2id 



Q/0 e 210 ,v l , ,„e l9 ,v l . ,„,v l „,„e' ie ] T . Here u and v 

n, — 6/2 1 n, — 1/2 1 n,l/2 1 n.6/2 ' 

are functions of r only. The special form of $^ (r) is cho- 
sen to preserve the particle-hole symmetry at I = and 
to remove the 9 dependence in the BdG equation ©. 
If $„(r) is a solution with an energy E, then there is 
another solution with the energy — E in the —I channel. 
Henceforth we only consider E > solutions. 

Generally the BdG equation © with a vortex can- 
not be solved analytically. Here we numerically solve the 
Eq. ^ and calculate the quasiparticle eigenenergies and 
eigenwavefunctions. In the calculation, we use the pair- 
ing gap A s from a self-consistence solution of the BdG 
equation for a pure s-wave superconductor with a small 
size of the system R = 25fc~ 1 (the Fermi vector fc c for the 
s-wave superconductor is chosen as 0.5 nm _1 ). Because 
the pairing gap approaches the bulk value in a distance 
much larger than fc" 1 , we can extend the pairing gap to 
a larger system R — SOOfc" 1 by inserting the uniform 
bulk value (see the inset in Fig. [3]). We find that there 
exists a unique zero energy solution when the chemical 
potential \i lies in the gap between the lowest two HH 
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Figure 4. Plot of the quasiparticle energies in a vortex core 
with respect to the chemical potential /i. Solid: zero energy 
state; Dashed: the minigap. Dotted: the bulk excitation gap. 
The parameters are the same as that in Fig. 1. 



bands (Fig. QJ>). In Fig. [3J we plot the two components 
Uq 1 / 2 (V) and WQ 1 / 2 ( r ) °f ^ ne zero energy wavefunction 



$o(r) and find w° 1/2 (r 



x 0.m 



(r) 



J 0,m. 



(r). We also confirm 
, (r) for other mj. Therefore the 



J 0,l/2 



Bogoliubov quasiparticle operator, defined as 

7n = * / dr Y] Kmj( r ) c L 7 ( r ) + Vnmj(r)Cmj (r)] , 

(10) 

satisfies 7q =70, which is a self-Hermitian Majorana op- 
erator. Consider two Majorana operators 7^4 and 73 in 
two vortices. It is easy to show 7^ — > 73, 75 — > —7^4 
upon an exchange of two vortices [20J . Therefore the Ma- 
jorana zero energy modes satisfy the same non-Abelian 
braiding statistics as that in a chiral p-wave supercon- 
ductor/superfluid [2l|, |22j and can be used for TQC. 

In Fig. |H we plot the zero energy state (the lowest en- 
ergy level at the I = channel), the bulk excitation gap 
(the first excitation at the I = channel), and the mini- 
gap energy (the lowest energy level in the vortex core at 
the I = 1 channel) with respect to the chemical potential 
fx. When \x lies in the gap between the lowest two HH 
bands, there exists a unique zero energy solution, which 
originates from the broken time reversal symmetry of the 
chiral f + if -wave superconducting pairing. The minigap 
is the topological gap protecting the Majorana fermions 
at the zero energy states and the associated non-Abelian 
braiding statistics from finite temperature effects. The 
numerical results show that the magnitude of the mini- 
gap is at the order between A s and A 2 /Ep- Additional 
numerical calculation shows that Majorana zero energy 
solutions also exist for a vortex with a winding number 
—3 or other odd numbers. We also find that the Ma- 
jorana zero energy modes exist when the second lowest 
band is LH, instead of HH. 

The existence of the Majorana zero energy modes can 
also be demonstrated analytically in the strong confine- 
ment limit. In this limit, the single particle Hamilto- 
nian Hq is replaced with Hhh in ([7]). The spinor wave- 
function at the angular momentum I channel changes 
to ¥ n (r) 



e ue \u l nt e- 



i0 n ,l 



J 



-lid 



The 



corresponding BdG equation can be further reduced to a 



2x2 matrix form 

Fo-M + 3/io I3Fi + \A s 
-{3F 2 -\A S Fi-fi-3ho 



Uf (r) 
H ( r ) 



(11) 



for a zero energy state after 6 is eliminated using 
the wavefunction & n (r) and the particle- hole sym- 
metry of the wavefunction is taken into account 

0. Here To W = Vo (Q ~ r~ 2 ) 2 + m (Q ~ r~ 2 ) , 
Fi(r) = d r (d r + r- 1 )(d r + 2r- 1 ), F 2 (r) = 
{dr + r-^dridr-r- 1 ), F 4 (r) = Vo(Q-4r- 2 ) 2 + 

Vi (Q - 4r~ 2 ), Q = ^ + ^ rl u a (r) = \ V<T (r). We ap- 
proximate the radial dependence of A s as a step function, 

1. e., A s — for < r < £, and Ao for r > £. Detailed 
analysis of the wavefunction (uf (r) ,u\. (r)) T shows that 
there are four and five independent solutions of Eq. 
(fTTj) inside and outside the vortex core respectively in 
the parameter region A = — 1 and p, 2 + A 2 < 9h^. The 
corresponding 9 unknown superposition coefficients for 
the total wavefunction match with the 9 constraints 
from the continuity of the wavefunction (up to the 
third order derivative) and the normalization condition, 
yielding a unique zero energy solution. 

In summary, we show that a chiral / + if -wave super- 
conducting pairing and the associate Majorana physics 
may be induced in a hole-doped semiconductor thin film 
through the proximity contact with an s-wave supercon- 
ductor. The proposed Berry phase mechanism in this 
system presents a new possibility for studying uncon- 
ventional pairing symmetry, which is distinctly different 
from the conventional scenario in which the pairing is 
induced by the Boson-exchange electron-electron inter- 
action mechanism. 
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